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Abstract. A quantum stochastic model for an open dynamical system (quan- 
tum receiver) and output multi-channel of observation with an additive nonva- 
cuum quantum noise is given. A Master equation for the moment generating 
operator of the corresponding instrument is derived and quantum stochastic 
filtering equations both for the Heisenberg operators and the reduced density 
matrix of the system under the nondemolition observation are found. Thus 
the dynamical problem of quantum filtering is generalized for a noncommu- 
tative output process, and a quantum stochastic model and optimal filtering 
equation for the dynamical estimation of an input Markovian process is found. 
The results are illustrated on an example of optimal estimation of an input 
Gaussian diffusion signal, an unknown gravitational force say in a quantum 
optical or Weber's antenna for detection and filtering a gravitational waves. 



Introduction. The time evolution of quantum system under a continuous ob- 
servation can be obtained in the frame work of quantum stochastic (QS) calculus 
of output nondemolition processes, firstly introduced in |5j and recently developed 
in a quite general form in [3| Q| Q] [5] . A stochastic posterior Schrodinger wave 
equation for an observed spinless particle derived in ^ by using the quantum fil- 
tering method .5 1 , provided an explanation of the quantum Zeno paradox [7| |5| . 
In this paper we give a derivation of the reduced wave equation for a Markovian 
open system described by Heisenberg quantum stochastic operators X (t) with re- 
spect to noncommuting Bose output helds Y (s) , s G M + which are assumed to 
be nondemolition in the sense 0J |3] of the commutativity [X(t), Y(s)] = at 
each time t > s. We shall obtain it by a non-unitary dilation of the characteristic 
operator of an instrument for the observable output process, but in contrast to [H] 
we restrict ourselves to the diffusion observation, i.e. to a continuous nondemo- 
lition measurement of a quantum Brownian motion. This gives the possibility to 
solve the dynamical problems of quantum detection and estimation theory |12| as 
demonstrated in an example. 



1. The dynamical model 

We are going to describe a dynamical model for continuous in time indirect 
nondemolition observation of an arbitrary family Q = (Qi, . . . ,Q n ) of Hermitian 

Date: September 20, 2005. 

Key words and phrases. Quantum input-output processes, Quantum nonlinear filtering, Quan- 
tum stochastic equations, Quantum optimal estimation. 

Author thanks to Matt James for support and hospitality in ANU, Canberra, where this old 
paper was revised and prepared for the current publication . 

This paper was read in the 1st QCMC conference in Paris, November 1990, and was originally 
published in the conference Proceedings: Quantum Aspects of Optical Communication, Lecture 
notes in Physics 45, pp 131-145, Springer, Berlin 1991. 



1 



2 



V.P.BELAVKIN 



operators Qj = Qj , acting in initial Hilbcrt space Ho of an open quantum system 
(antenna), with additive ^-correlated Bose type quantum error-noises in linear n- 
dimcnsional output channel of observation ej(t), j = 1, . . . ,n. We shall describe the 
output quantum error-noise e(i) = (ei, . . . , e„) (i) by the components Vj (t) = ej (i) 
of a quantum stochastic process which is opposite, or inverse to an input quantum 
noise v. = (vj) of the same or higher dimensionality m > n. The total quantum out- 
put noise v. = (vj) as the opposite to v. can be defined as a compatible (commuting 
with v.) but maximaly closed (maximally entangled with v.) quantum stochastic 
process v. which simply coincides with the input v. if it is is self-compatible, i.e. is 
classical (having all commutig components Vj). 

Let us describe the quantum input noise v. by the classical white noise compo- 
nents Vj(i) represented by noncommuting Hermitian operator- valued distributions 
Vj = Vjt. They are completely determined in a quantum Gaussian state by the first 
and second moments 

(1.1) <Vi(t))=0, <Vi(t)v fc (t')> =K ik 6(t'-t). 

Here Ki k are complex elements of a Hermitian-positve matrix k — [nik] of the 
same or higher dimensionality m > n, with imaginary part Im k defining the Bose 
commutation relations 

[vi (t) , v fe (t')] = 2ilm K ik S (t - t') 1 , 2iImK jfc = n ik -K ik , 

such that complex conjugate components define the intensity covariance 

matrix k — [nji] of the output noise v. = (vj) : 

(vi(i)vfe(i')) = K lk S(t - t'), [vi (t') , v k (t)] = 2iImK ifc /<5 (t - t') 1 

as transposed (or complex conjugate, 7c = k) to k. Thus all output components Vj 
commute with all input components vfc, and v must also be maximally correlated 
with v in the sense that the intensities jj k of real covariances 

(vi(t)v k (t')) - lik 5{t t') = lki S(t' -t) = (v fe (t) Vi (t')) 
are the elements of a symmetric m x m-matrix 7 = [7^] as the geometric mean 

~ 1 /2 ~ 

7 = (/t • k) 1 . The geometric mean with k for an invcrtible k is defined as a 
Hermitian-positive matrix 7 such that k — jk^ 1 ^. The matix 7 is symmetric and 
invertible, with the inverse 7 _1 defining kT 1 = 7~ 1 K7~ 1 and as the intensity 

for the covariances 

(v l (t)v k (t')) = n lk 8{t - t'), [v l (t) , v k (t')} = 2iIm K lk 5 (t - t') 1 

of the contravariant components v l (t) = 7 v v. (t) for the output noise v. (t). (We 
assume that Hermitian matrix k = [nik] is strictly positive, with the inverse k^ 1 — 
[n lk ] corresponding to a finite temperature of the output quantum noise v' = (v 1 )). 

As usual the operator-valued distributions v, (t) can be described as general- 
ized derivatives Vj (t) — ^v* of quantum Wiener vector-process vf represented by 
selfadjoint quantum stochastic integrators v* = v'^ with = and independent 
increments dv* = v' +dt — v' satisfying the noncommutative multiplication table 
dv*dv| = n ik At. Assuming for simplicity that matrix n commutes with the trans- 
posed k one can realize such quantum Wiener noise with respect to the faithful 
state given by the vacuum vector 6$ in a Fock space J- as v' = a* +&y = 25Ra* , the 
doubled Hermitian parts of the linear combinations a* = K 1 / 2 a t . Here & k are the 
canonical annihilation integrators which are adjoint to the creation operators a'-* 
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in the intervals [0, t) defined on the symmetrical tensors over the complex vector- 
functions a.(t) — [at, Oil-, ■ ■ ■] (t) with 

(a. | a.) = YL J "(*) a (*) d * =11 a ll 2 < 00 

■ik 

as the symmetric tensor multiplication by the indicator function l[o.t)- The canon- 
ical commutation relations 

(1.2) [a(a*), aV)] = H«0 . [a(a.)> a ( a ')] = 

then are realized by the quantum stochastic integrals a(a*) = J afc(t)da t , a? (a.) — 
&(a*y . (We use Einstein notations for the convolution ak/3 k = J2&k(3 k over the 
indices k — 1,2,... in contrast to the scalar product notations /3 • a* for the finite 
sums Y^j=i P 3 (Xj) an d omit the identity operator 1). 

The output vector-process Y(t) = Q(t) + Jo (8 e(t), defined by the integrals 

(1.3) Y j (t)= [ Qjtfdr + IoQel, j = l,...,n 

Jo 

of the Heisenberg operators Qj(t) = U(ty (Qj (g> 1) U(t), can be realized for a sin- 
gular coupling of the system with the Bose fields a| by the output observables 
e* = 23?aj, a* = K 1 ^ 2 a t in the interaction picture 

(1.4) Yj(t) = U{tf (j„ (» e*) U(t) = 2$tBj(t), 

where Bj(t) = U[ty{I <g> a*-)J7(i) are the annihilation output processes, introduced 
in ^1 and Jo is the identity operator in Ho- The unitary evolution U(t) 

will be described on the tensor product TL = "Hq®3 : k , by a Schrodinger-Ito quantum 
stochastic equation 14 

(1.5) dU{t) + KU{t)dt = i Q Q ® d/ (i? t ) - 29(4 (8 da£)^ 17(f), 

in terms of the input integrators h t — k _ — 7~ 1 a t , where Lk, fc = 1, 2, ... are 
the operators in TCq with Lj + iA = Qj, 7 = 1,..., n, 

/'($) = 4^ /($), -ff = is a Hamiltonian of the system, and Q — is an operator 
in 7Y of a generalized coordinate conjugate to the generalized force / (f) = ^/ (i? t ) 
depending on an independent input diffusive signal $ t , the random position of a 
gravitational source say, with (d$ 4 ) 2 = <r 2 dt. Note that in the case Lk = L\ this 
equation can be written as 

dU(t) + KU(t)dt = ~ ( Q ® df (0 t ) + Q fe ® df t fe ) l/(t) 

in terms of the integrators f* = — ft^a^ of quantum Langevin forces f fc (i) = ^f t fc 
satisfying the canonical commutation relations 

[P'(t), ej -(t0] =y^(*-0. FW, f fc (*')] = if (*-*') 

corresponding to noncommutative multiplication tables 

dejdf? = ihS k di, df 'df? - (ft/2) 2 K lfe dt, 
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including S*- = {2a f /h) 2 S l Q for f° = / ($ t ). 

The solution t i— > U(t) of the equation (|1.4(l for t/(0) = I is adaptive C/(s) = 
U S ®I[ S with respect to the tensor decomposition TL = Hs^J 7 ^, where 7i s = TLq^Ts 
and !F a , are the Fock spaces, generated by the vector-function with a s (t) = 
0, Vt > s and at[ s (t) = 0, t < s respectively. Hence (I <S>ej)U(t) = U(t)Yj(s) for 
any s <t, and 

[Yj(s), Y k (t)} = U(tf(I ®[2m%, 2$t&i])U{t) =2isImK jk I 

(1.6) [Yj(s), X(t)} = U(t) f [I Q ®ej, X® l]U(t) = 0, Vs < t 

for any operator X(t) = U(ty{X ® l)U(t) of the system in the Heisenberg picture. 
This proves the nondemolition property of any observable process Yj (t) with respect 
to the system. It was introduced as the nondemolition causality principle for the 
output quantum processes in E] • Using the quantum Ito formula (| 1 . 41) for 4 
with dej = 2Kda,j and multiplication table 

(1.7) dajdaf = n tk dt, dajda^ = S^dt, da^da^ = K ik dt 

with all other products being zero, one can derive the equation l|1.3[) : dY(t) = 
Q(t)dt + J <g> de(t), where Q(t) = [Q x , . . . , Q n ](t), Q 3 -(t) = U(t)t(L 3 - + L))U{t). 

Theorem 1. Let us suppose that the input real-valued signal $ t satisfies the sto- 
chastic differential equation 

(1.8) d&t + v(& t )dt = adw t , 

where w t is an independent standard Wiener process, defined by the moments: 
(w t ) = 0, (w s w t ) = s for any s < t with respect to the vacuum state as w t = 25fta° 
for the canonical annihilation integrator a° in the Fock space T . Then any twice 
differ entiable in the strong operator topology function I:i)E!Rh X({t) G B{TLo) 
in the Heisenberg picture X(t) — X*(i9 i ), where X 1 — U (t) XU (t), satisfies the 
following quantum stochastic Langevin equation 

dX(t) + (vSX{t) + ±[X{t),H (t)] - X - (a 2 S 2 X{t) + A? [X (t)])^ dt 

(1.9) = [X(t), 23?4] <g> ^df t fe + [JSf(t), iQL\] ® dv t fe + (j SX (£) ® dw t . 

/fere vj = 2Kaj = 7 * fe v£, <SX(tf) = X'(tf) + /' (0) [X(<&), j-Q], X'(0) = £ X(0) 
and 

(1.10) Af [X^)j = i f] [X t W , 4] + [ i *t ) jf'(0)]L* fc ). 

i,fc>l 

2. The reduced evolution 

Let A denote the input-system algebra which is assumed to be the von Neumann 
algebra of all essentially bounded operator- valued functions X : i3 e R i— > X($) £ 
B(Ho), bt C B(Tt) be von Neumann subalgebra generated by the error-noises 
{ef,...,e*} for all s G [0,t), and <? t C :F f be subspace generated by b t on the 
vacuum ^0 6 .Ft for each i > 0, where Tt are Fock subspaces generated on 8® 
by all input processes v*, w s or, equivalently, by forces (fj,f *,.••) up to time 

t. The increasing family |bj : t > o| with respect to the output states induced 
by ip(t) = U (t) (tp ® 5$) S Tt, -00 e ^0 will be called the output filtration in 
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the Schrodinger picture; It is equivalent to the Heiscnberg filtration {Bt,t>0} 
with Bt C A ® B (Ft) generated by the output family {Y\ (s) , . . . , Y n (s)} for all 
s E [0,t) with respect to the initial states induced by ip (0) = Vo ® ^0> smce 
[7 t 'l^- (s) U t — with {J { ' = [/ (t) |Wt for any s < t. The filtered dynamics of a 

quantum stochastic system, described in the Heisenberg picture by homomorphisms 
X (t) = U t XU} of A 3 X into A ® B (F t ), is defined by the cocycle of CP maps 
$i : .4. — * A® It, where l t are (unbounded) commutants of b t on Q t , such that these 
dynamics induce the same input-output states on A® bt with respect to the initial 
vacuum state: 



(2.1) 



C0 



[*] (*> 



® b 



f 



J7 t X(tf t )®6 C/ ( 



for all X e 4. and 6 E b t . Here eg [•] = (Iq ® 8®) [•] (Iq ® 8®) is the vacuum (condi- 
tional) expectation A® B (J-) — > A such that the composition p o e® = e® o p with 
any normal state p on A is the product state p ® e® . 

Since each $ t is normalized as <& t (la) = Pt to a positive element P t G ^4® 
defining typically unbounded density operator pt — p [Pt] £ £ for the output state 



U t (lo (8 6) C// 



£0 



bp ( 



ft (b) = P (e® 

with respect to the input state e® b on the algebra bt, it is necessary to give a more 

precise meaning of the unbounded commutant £ t of b t . Due to the nondegeneracy 
of the covariance matrix-function [Hjj> min {s, s'}] of e s : s < t for each t > 0, 
the cyclic representation bt\Gt on the minimal invariant subspace Q t containing 
8® is faithful in the sense that b — in b t if b8® = in Ft- Therefore b t is 
transposed to its bounded commutant b' t , coinciding with bt = JbtJ, where J is an 
isometric involution defining the transposition bt — > bt by b' t — Jb\j common for 
all subspaces Qt- Moreover, the von Neumann algebras bt and b' t are in one-to-one 
correspondence with the achieved Tomita algebras bt = bt8®, b' t — btS® as 
dense subspaces of b — b 8®, b' — bS® in Q t , where b — JwJ, b £ b t , with common 



identity 1 := 8® 
II b |U:= ll&ll = 



involutions b" := wS®, b' b := b^Sa such that b"' = b'^ and the norm 



V ||c 



Let us define a dual space It to the Banach algebra b' t as the completion of 
bt C 5t with respect to the dual norm 

|| a |U= s U p{|<6', a ) | : & £ b' t , || b' [[oo<l}=NL. 

where the bilinear form is the standard pairing of bt and b' t , 



(2.2) 



(b',a) := e® 



hit 



b, a 



which we will extend by continuity on all a G It 2 bf. Thus the co-algebra £ t of 
the operator algebra b t is defined as the Banach space of operators a : b' i— > ba, 
mapping b' t C Q t into £ t D Q t , bounded with respect to the dual norm || • || ^ . This 

space is a linear span of positive cone It — {P> 0} defined by (V f b,p^ > for 

b^pb of b t 



all b > 0, and therefore is invariant under the right and left action p 
which is defined as the dual to the selfaction q\-*b'qb on b t , 



(2.3) 



(q,tfpb) =(tiqb,p) Vqeb t ,peb 
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extending the transposed selfaction on bt by b — Jw J for all b £ bt . The coalgebra 
it is also equipped with involution a >—> of defined by a*b' — ba^,Va £ it such that 
(b', a") = (b, a') id, and with the identity 1 = 1, corresponding to the vacuum state 
e = -Dfl- Note that the elements ii6< obviously commute with b € bV 

a&c' = 6ca = feac' Va £ £t,b, c £ bt, 

but are unbounded in the Hilbert space Gt if o ^ bt. However they are densely 
defined as the bounded kernels in the Gelfand triple bt C (J t C £ t , and = 
(l,p)g for any positive element p € i t which means that it is density operator of a 
normal state on b t if (l,p)^ = 1. Moreover, every normal state is uniquely given 
by such density, i.e. that the space it is predual to bt, and is preadjoint to bt 
which we denote as t\ = bt, i\ = bt- In most cases the density operator p £ i 
of an output state q = (j>,P^ has the range pi in J- t , as it is in the case 

of a majorized positive form s(wb) < Xe^wb) for a A > 0, corresponding to the 
bounded p on J- t : \\p\\ < A; moreover, any operator p £ i can be approximated by 
the density operators from the bounded commutant bt in the norm =|| p ||i, 
where p — pSy, . 

In order to derive a quantum stochastic equation for the reduced dynamics 4> t , 
let us find the differential evolution for the factorial generating map : A — > A, 



(2.4) =e [$t[X](Z®^)] = (8 t \xW\t)\5 9 ). 

Here &t = zf^ are exponential elements, defined by the Wick (normal) exponent 

(2.5) zf ] = t&WJfMPt) = : exp [ et (/3)] : 

of the observable e(/? t ) = a (/3 t )+a^ (j(3 t ) = Jo /^( r )'de r , where /3 = is a column 
of locally square-integrable functions with (3 3 8 {t) = for t > s, j > n, (3 J s (r) = f3 3 (r) 
for r < s, and X^(t) = U t {X{d t ) ® z[ P) )ul = X(t)Z(t). Taking into account the 
equation l|1.9|l and 

dZ(t) = Z(t)/3 j (t)(I ® 2d5Ra* + U t QjU}dt) 
for Z(t) — Ut (io <S> if' 3 ' 1 ^ C^t , one can obtain by quantum Ito formula ^1] 

d(X(t)Z(t)) = U t (A*[X}(#t) + (tfX^ + X^L)P(t))u}dt 

+aU t SX^ ) u} ® dw t + U t X^U}0 j (t) <g> d25fta* 
+£/*[££, ® da t fe + da^ ® C/ t [Jff } , 

where AT t (/3) = X(0 t ) ® z t (/3) , 6 = 8/&& t + f [; {Q\- Hence the map $f ) : X £ 
A i ► ((50|J7fXj /3 ^[//|(50) satisfies the equation 

(2.6) ±$\ f}) [X]=$[ l3) [A*[X] + (tfX + XL)f3(t)}, ^ [X] = X , 
where A* [X] (0) = -v (0) SX (0) - | [X (0) , if] + A (a 2 5 2 AT (0) + AJ [X (0)]) . 
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In the same way, using the Ito formula for the product GtZt of Gt — $t[-X] and 
the Wick exponent Q2.5JI . one can obtain the equation for (|2.4(l if 

n 

(2.7) d$ t [X] - $ t [A*[X]]dt = ML]X + XLj] ® dv?, 

where the operators v| = 25fta,(, satisfying the canonical commutation relations 
[vj , Vj,] = 2iIniK ; ' J min{£, i }, 

generate the predual coalgebra^t as the unbounded commutant of bt = {ej : s < t}" 
since v^,j = 1, . . . , n leave the subspaces Qt invariant for any t > s, and {ej : s < 
iji < = bt = {v^ : s < t, j < n}" on Q t . This proves the following theorem. 

Theorem 2. Let the initial state p : A —> C be a normal one, described by a 
density g : 1? i— » f?(#) wii/i values in trace-class operators on Ti.o such thatp (X) = 
J r Trg('3)X (??) d&. Then the conditional state 

(2.8) h[X\ = J Tr H M^t[XmW 

is described in the standard representation by the operator-function (p t ($) as <p t [X] = 
fTr-n [ip t ('&)X('d)]d'd satisfying the quantum stochastic equation 

n 

(2.9) dft(0) = A[ft](0)dt + £>^ t (i?) + ® dvf 

wit/i <^o(^) = -ffere the quantum stochastic differentials dv J t together with de* 

have the canonical multiplication table 

dv{dv{ = K jj 'dt, dv^'de*., = 5 j f dt, de*dej, = «#/di, 

and A [99] = 5 (vip) + 4 [<p, if] + | (cr 2 (J 2 !/? + Ai [y>]) is preadjoint generator defined 
on A (1?) - ^ (1?) + /' (tf) [v, (1?) , 1Q] , 

<5 2 ^) = + (2f (0) +f"($))&i(d),^Q}+f(d) 2 [W)^Q],^Q], 

(2.10) Ai [tf (0) = J2 * ik ([Li, + [J^W, 41). 

i,fc>0 

If is normalized to a positive martingale pt — J Tr-}{ ip t ('&)d& £ £ t as i/ie density 
operator defining by \2.1\) the output state (b(t)) = p[(5$\B (t) \6q)] on £' 3 z for 
B(t) = U(ty(I ®b)U(t) with respect to the vacuum state vector 5% £ T simply as 
(S(t)) = (ft,6) . 

Note that since all the input components v\,i = l,...,n commute with the 
output components e*- and have with them zero correlations and thus are indepen- 
dent of unless i = j, they generate the same subspaces Qt as e*-,j = 1, . ..,n. 
On these subspaces they simply coincide with the transposed ef = Je l t J = to 
the contravariant components e\ = 25RaJ of the output noises et, having the same 
autocorrelation and mutual correlation matrices with e k as vj and being given by 
k t = K _1 / 2 a t = 7 _1 a* such that e\ = Y^k=\ l %k e k ■ Thus in the filtering equation 
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(|2.9(l the input noises n\ on Qt can be replaced by e J t ' which can be given in the 
Gaussian case as linear combinations of the transposed noises e'j — Jej J = e'j . 

3. The optimal measurement. 

The quantum filter defines a quantum measurement on the output of the system 
at a time instant t £ R + . In general it is described by a i?t-valued positive measure 
M(t,dx) on a Borel space X, normalized to the identity operator in H = Ho <8> 
T: J M(t, dec) = I. The problem of optimal quantum observation is the problem of 
finding the optimal measurement M°(t), giving the minimal value of the mean 

(3.1) J (S(t,x)M(t,dx)) = J (m' t (dx),c t (x)) = J (m t {dx),c t {x))$ 

for an „4-valued measurable function S : x e X i— » S(x) — S(xy in the Heisenberg 
picture S(t, x) — UtS(x, $t)U} with respect to an initial state (•) = p o e$ [•] on 
A <g> B t - Here the mean (|3.1fl is given in the standard representation B t — > bt 
according to Ij2.1|l as the integral of the pairing (|2.2|l for 

(3.2) rn' t (dx) = m t (dx)5y > , c t (x) = c t {x)Stj S , 

where rht defines the measure M{t) — Ut(I <£> rht)U} in the Schrodinger picture, 
and 



(3.3) ct(x) = J Tr Ho [ip t (#)S(x,#)}dti = <t> t [S t (x)]. 

The duality principle gives the necessary and sufficient conditions |13| of optimality 

(3.4) J (m t (dx),c t (x)} = hff { J{m t (dx),c t (x)) : j m° t [dx) = 1} 

of a positive bt-valued measure vn\ > with J rh°(dx) — 1 as the conditions for 
the dual problem 

sup{(ij) : l t < c t (x),x 6 X] = (iJthr 

where It — 1° defined by kb' — bl t is the operator 1° < c t (x), Mx E X, for which 
/ (rht(dx), Ct{x))% = (1, The last condition of optimality can be written in the 
form of the equation 

(ct(x) - lt)m't(dx), or m t {dx){c t {x) - l t ) = 0. 

Let us consider now the problem of optimal estimation \T2\ of a selfadjoint op- 
erator X(t) = JT'(t9 t ), given in the Schrodinger picture by a B(TCo)-valued function 
X(<&) = X(-i?)t, with x = A e R, S(x,0) = (X(ti) - xlf. One can treat in such 
a way the problem of filtering of a real measurable function x("d t ) of the input 
diffusion signal ■dt taking X(d) = x($)I. In order to formulate the optimal esti- 
mate in terms of the measurement of the optimal output observable x^ € bt as an 
appropriate posterior mean of X (t) we need the quantum stochastic equation for 
4> t [X] = J(j) [X'\ J in terms of the output noise e*. It can be obtained by complex 
conjugation 



n 

(3.5) d& [X] = ^oA* [X] {$)dt + J2^t [LjX + XL] 



>de 3 t 
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of the equation for the CP map <p t : A — > bt described in (|2.8(l by the density 

operator <p t satisfying the equation i|2.9|) . Here ej. — 53j=i ^ e j'i where 9" 

9 1 is real symmetric n x n- matrix with the inverse 9 = [Ojj>] as an intensity matrix 
of the standard covariaces 

MOe^t)) = e rj 5{t' -t) = e J3 ,5{t - if) = ( c ;- (t) e f (0) 

of the output noises &y with the transposed components e'j = Jej J = e'J such that 
9k l 9 T = k for the n x rt-submatrix k of k. (It is the geometric mean of k and k, 
e.g. 9 — k^7«2 if k and 7t commute.) Note that since in general 9 is smaller than 
the n x n-submatrix 7 = [7,y] of the mutual covariance matrix 7 for e.; = Vj and 
Vfe, e^- 7^ Vj unless m = n, and similar e J t ^ v J t since ej = v t on and 7 _1 < 

Theorem 3. The solution of the optimization problem for the quadratic cost 
function 

c t (X) = X 2 p t -2X^ t [X} + MX 2 } 

is given by the spectral measure fh° t of a selfadjoint operator x\ E bt defined by 
— Jx t J as transposed to the operator xt = x\ resolving in bt the equation 

(3.6) x t p t +p t x t =2 J Tr-tolX 

It is given as the symmetric posterior expectation x° — p t [X] by an operator- 
function 1? 1— > g t (1?) E B* (7i) (g> It with It = J it J as the density of the solution to 
the equation 

(3.7) ptpt [X] + Q t [X] p t = 2$ t [X] VX E A, 

where pt = JptJ ■ The symmetric posterior density g t = g[ satisfies the nonlinear 
quantum stochastic equation 

n 

(3.8) dM#)=rt[Q t ](m + Y, E i^® de % &dW = eW- 

where T t (t9) and 3j(l?) are defined by the solutions to the equations 
UP0[g}{'&) = SR (Pt&QW + - P 3 tQ{$)) , 

n 

(3.9) mmm = u(p t A[Q](#)-Y, Kjip i E M)> 

with Jj = Y%=i L ^ ij ' Pt = h® Pt and PI = I ® pj defined by Q j = V + and 
Qt as 

(3.10) f t = Klptqf], 4 = J Tr Ho [Q t WQ j W- 

Proof. Denoting u t — J Arht(dA) for an orthogonal projective- valued measure 
TOt(dA) E bt, one obtains 

<m t (dA),c t (A)} = (ulpt)-2(ut^ t [X]) + (lMX 2 ]) 

= (ul,p t ) -2(K(x t Ut) ,Pt) + (l,MX 2 }) >(U?>, 
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where it is defined by the duality (|2.3J) as the transposed to the solution x of the 
equation l|3.6|l written as 5Rp t it = </> t [^], and 1% £ it is defined as 

l° t = fa[X 2 ] - itptXt- 

The inequality is due to positivity of pt — 4> [Iq] £ it and (u — i) 2 £ b*: 

<m t (dA), c t (A)) - (i, 1°) = ((u t - x t )\ p t ) > 0, 

and the equality is achieved at the spectral measure m° t of the selfadjoint operator 
it £ bt defining u t as it- 

Representing the solution of the equation Q3.fi|l as 

x t = J Tr Ho X($) g t (#)d$ = p t [X] 

in terms of the density function g t £ A* ® it we note that p t satisfies the equation 
transposed to (|3.7|) . and therefore p t [X] = Jp t [X^] J satisfies the equation (|3.7|> . 

Looking for the operator i t — p t [X] as the solution of a quantum stochastic 
equation 



di t = g t dt + c{ ® de 



j 



with some gt = g\, & t ' = H > we should compare the quantum stochastic differential 
for </> t [X] = J~<t> t [X^] J, with the differential 

n 

dStptXt = U(dp t dx t + dptxt + ptdx t ) = Kji(f> t [Q l ]c? t dt 

n n 

+5R MQ l \xt O de- + 5Rp t (^ c> t ® de* + g t dt), 

i=l 3=1 

for tftpt&t = <l>t\X\ obtained applying the quantum Ito formula, where we took 

n n 

dp t = k\Q' ] ] ® de*-, Q j = Qift 

j = l i=i 

for the martingale p t — <j> t [Iq]- Here e| = Y%=i ®ij e t as we have expressed the 
differential dpt — Y^j=i fa[Qj] ® ^ e t f° r the martingale pt = Jpt J in terms of the 
driving output error noises e*- by the real linear transformation 6~ l . Comparing 
dSiptXt with the equation (|3.5|l for d(j> t [X] 



[X] = fao A*[X](0)dt + ^ fa [ vx + XLjf ] ® de 



/ 

j 



written in terms of e*- , we derive 

M[fa[Q3]x +p t c k t ] = fa[L^X + XV] = ZtptpAV^X + XV], 
^faun+ptgt] = fa[A*[X]] = np t p t [A*[X]], 
where p t [X] = Xf This gives (|3.8|) - (|3.1U|) in terms of 2j, T t , defining 
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Note that the unnormalized posterior expectation <j> t : X i— » <f> t [X], as well as its 
normalized version p t [X] — jTt-H [X (fl) g t (d)]dd is not CP map but transpose- 
CP, in the contrast to the CP map (f> t : A — ► bt and its normalized version p t [X] = 
Jp [Jft] J. 

Example. Let us consider the case n = 1 with L\ = h Q = L\, 7 1 ^ 2 at = a\ — 
7~ 1 / 2 a t , where a.] is the standard annihilation integrator [aj, &l ] = min{s, t} defin- 
ing the input noise v t = 25ftat and nondemolition observation of the commutative 
output process 

dY(t) = Q(t)dt + I ® de*, Q(t) = U{t)^QU{t) 

by e* = 23?a* = 7v t . In this case et = v t , y* := e'J© is v' t — jvt, where vt — vtS® 
is given as j~ 1 ^ 2 wj by the standard Wiener process w\ identified with the vector 
process &\ 5$ in Fock space, and the output process Y (t) on the initial state vector 
t/>o<8><$0 is identified with the classical output process y l = rf^w] = jvt relatively to 
the probability density p t (v) = ptS® with respect to the Wiener probability measure 
P 7 of the input Wiener process {v t } with the intensity 7 . The equations (|3.5|) 
and l|3.8|) are classical stochastic equations in the linear (for the nonnormalized 
ip t (i!},v) = <^ t (z?)^0) and nonlinear (for g t (i!},v) — tp t ( j d,v)/p t (v)) posterior density 
operators with respect to the output states 

(V»(8 h\B t ip® 6 9 ) = / b t (v)p t («)P 7 (d«) = (6t(v),Pt(v)) 



on B (t) = b t (7 Y) given by any adapted functional bt of v. 

Suppose that the quantum receiver is an open oscillator (e.g. Weber's antenna), 
described by the Hamiltonian H = h A^A, where A = iP + ujQ and Q,P are the 
canonical coordinate and momentum operators: [Q, P] — ihl . Then the quantum 
Langevin equation 11.9(1 for A(t) = UtAJj} is the linear one 

dA{t) + iuA(t)dt = U <8> {dd t + df t ), 

where f t = SSaj is the Langevin force (thermal noise) as a classical Wiener process 
of the intensity er 2 = tf/ij acting on the coordinate Q, defining the total force in 
the right hand side of the equation as the sum $ t + f t of the unknown gravitational 
force / ($t) — &t an d the thermal noise through the additive channel d t (w) + f t . 
In the case of Gaussian input process i?t, corresponding to the linear v(i9) = v$ 
and Gaussian initial state <f> the optimal estimate of X(-dt) = $tl = X (t) is given 
by the linear posterior mean value $ t = it with respect to the output coordinate 
process Y(t). In the standard Fock representation it is given as the last component 
of the stochastic row x t — [qt,Pt,&t) 01 the posterior mean values for X(t) = 
(Q(t), P(t), X(t)), satisfying the Kalman equation 

dx t + x t Adt = k t dv t , A = 

where k t = (kj 1 , fc t 12 , k} 3 ) is the first row of the symmetric 3 x 3-matrix K\ = (k\ 3 ) 
satisfying the Riccati equation 

d 1 / 

— K t + K t A + A T K t H — kjk t — T, T= a 2 + a 2 a * 

dt 7 




a 



2 „2 
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with an initial symmetric covariance matrix Kq of (Q,P,$o) and vt = vt — 7 
The pair (xt, Kt) defines the posterior (normalized) Gaussian state of the quantum 

system with input signal Xt = fit and the mean square error ($ t ) , i? t = § t — $t is 
given by the component kf of the posterior correlation matrix Kt ■ 

A posterior dynamics of the quantum system under another nondemolition mea- 
surement of the received electromagnetic field by a photon counter is considered in 

m 
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